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I Abstract 

^ ■ We formulate the perturbative QCD approach to coherent diffractive dijet produc- 

00 . tion in pion-nucleon and pion-nucleus colhsions at high energy. For hard dijets the 

^iJ \ Ponieron sphtting mechanism in which both hehcity amphtudes are proportional to the 

Qv^ ■ unintegrated gluon structure function of the proton J^{x, k ) and pion distribution am- 

^D . plitude (pniz) is shown to dominate. In nuclear diffraction multiple Pomeron splitting 

^^ I components are found to give antishadowing contributions at large jet momentum k. 

■^ ' To leading twist there is an exact cancelation of effects of nuclear attenuation and anti- 

f^ ■ shadowing/broadening of multiple Pomeron splitting contributions. The next-to- leading 

' i higher twist correction driven by nuclear rescatterings is calculable in hard QCD and 

g_) I proves to be numerically very large. We argue that large higher twist effects do not 

r^ • preclude the determination of gross features of (I)tt{z)- Our results on the atomic mass 

^ , number and momentum dependence of dijet cross sections agree well with the preliminary 

k>( I findings from the E791 experiment. 



1 Introduction 

Ever since the classic work in early 50 's by Landau, Pomeranchuk, Feinberg and Glauber on 
diffraction excitation of deuterons [|1], |^, ^ into the proton-neutron continuum the momentum 
spectrum of excitation products (protons & neutrons) is known to be given by the momentum 
distribution of constituents in the deuteron. More recent work on diffraction dissociation 
focused on diffractive deep inelastic scattering (DIS). Here the microscopic QCD description 
of diffractive scattering by exchange of a color-singlet two-gluon tower in the t-channel reveals 
a sensitivity of the mass spectrum in diffractive excitation of the continuum qq states to the 
gluon structure function of the target 0]. Furthermore, extending early considerations in [Q 
Nikolaev and Zakharov have shown in 1994 [^ that in diffraction excitation of hard dijets 
7*p -^ p'qq there exist two regimes depending on how the large transverse momentum k of the 
jets compares to the hard scale of DIS, i.e., whether k <Q or k >^Q. In the first regime the 
transverse momentum k of jets comes from the intrinsic momentum k of the quark (antiquark) 
in the qq Fock state of the 7* and diffractive amplitudes are proportional to the familiar 
integrated gluon structure function (GSF) of the target proton G{x,k'^), see also |^. In the 
second regime, k >^ Q, diffractive dijets are a unique probe of the differential (unintegrated) 
gluon structure function (DGSF) of the proton jF(x, Q^) = dG{x,Q'^)/d\ogQ^. Specifically, 
in this regime the transverse momentum k of jets is provided not by the momentum of q and q 
in the virtual photon, but by the momentum of gluons in the Pomeron. Correspondingly, this 
regime of diffractive DIS has been dubbed 'the splitting of Pomerons into dijets' 0. In this 
regime the diffractive dijet production amplitude is proportional to J-'{x,k'^). Subsequently 
Golec-Biernat, Kwiecinski and Martin |^ reformulated the formalism [^] in terms of the off- 
diagonal (skewed) parton distributions. Because the skewed distributions can be approximated 
[^ by the diagonal ones at a rescaled x, after this rescaling the formulas of |^ are recovered. 
In the present communication we extend the approach f^ to coherent diffraction of pions 
into dijets on the nucleon and nuclear targets. The principal novelty compared to photo- 
and electroproduction is that the pion-quark-antiquark vertex is non-pointlike which makes 
splitting of the Pomeron the ever more important mechanism for hard dijets. We focus on 
coherent diffraction production on nuclei which has recently been measured by the E791 
collaboration 0. The principal issues with the hard QCD interpretation of these data are 
whether 1.25 <k< 2.5 GeV is sufficiently hard for the pQCD treatment, how large are next- 
to-leading twist corrections, what are nuclear effects and whether the extraction of the pion 
distribution amplitude is possible from the E791 data. The two major nuclear effects one has to 
deal with are nuclear attenuation and nuclear broadening of jets. The practical calculation of 
diffraction on nuclear targets involves evaluation of multiple gluon exchanges between nucleon 
and excited qq system and we take full advantage of the recent determination of the DGSF of 
the proton |iy]. We demonstrate that the broadening of the jet momentum distribution comes 



entirely from the multiple Pomeron splitting diagrams. The large k behavior of T{x^ k ) found 
in 



p!0[ is shown to entail a remarkable cancelation of the attenuation and broadening effects 
to leading twist. In view of these cancelations the principal nuclear effect is a higher twist 
correction which is perturbatively calculable and is proportional to G'(x, fc^). 



The further presentation is organized as follows. In section 2 we introduce the principal 
formalism starting with excitation of diffractive dijets on free nucleons and isolate the two 
helicity components of the diffraction cross section. We demonstrate how the dominance 
of the Pomeron splitting mechanism into hard dijets and therefore the proportionality of 
diffraction amplitudes to jF(x, k^) do emerge because of the non-pointlike pion-quark-antiquark 
coupling. For the same reason diffractive amplitudes are shown to be proportional to the pion 
lightcone distribution amplitude 07r(-2) of much discussion in the recent literature (for the 



review see 



11. 12 



The possibility of measuring (f).n-{z) in diffraction of pions into hard 
dijets has been mentioned in |]l4l but as we show the claim in Jl^ that to the leading twist 
diffractive amplitudes are proportional to G{x,k ) is in error. Calculation of multiple gluon 
exchange in diffraction off nuclei to leading and higher twist is described in section 3. The 
novel feature of nuclear diffraction are multiple-Pomeron splitting processes in which the k 
distribution is broadened by the gluon momentum coming from different split Pomerons. In 
the standard nuclear multiple-scattering expansion the higher order nuclear rescatterings are 



known to generate nuclear shadowing |T^ . We demonstrate that after reformulation in terms of 
multiple Pomeron splitting components the nuclear multiple scattering expansion takes a form 
in which higher order Pomeron splitting components give antishadowing contributions, i.e. an 
enhancement of the corresponding impulse approximation term. We find an exact cancelation 
of effects of nuclear attenuation and antishadowing/broadening of multiple Pomeron splitting 
in the leading twist k distributions. The higher twist correction rises with the multiplicity of 
split Pomerons and is shown to be proportional to the integrated gluon structure function of 
the proton G{x, k^). It is an antishadowing correction and rises with the nuclear mass number. 
In Section 5 we summarize our main results and present a comparison with the preliminary 
experimental findings from E791. Our numerical analysis shows that the leading plus next- 
to-leading twist asymptopia sets in only for k > 2-3 GeV, somewhat beyond the kinematical 
range of E791 1.25 <k< 2.5 GeV. Our numerical results for the k and atomic mass number 
dependence of the dijet cross section are consistent with the experimental findings by E791 

i- 

2 Microscopic QCD mechanism of diffraction into dijets 



We only need a slight adaptation of the formalism developed in [|, |^, |T6|. Diffraction dissoci- 



ation of the pion into the high mass continuum, hard, qq dijet final state. 



J ^.^ 



TTp ^f p qq , 

is described by the four pQCD diagrams of Fig. 1. In this paper the dijet cross section is 
calculated at the parton level. The relevant kinematical variables are shown in Fig. 1, A is 
the transverse momentum of the excited dijet, quark and antiquark jets carry a fraction z and 
1 — z oi the pion's momentum and the invariant mass of the excited pair M is given by 



where ruf is the quark mass. Such a parton level modeling of final states is applicable if the 
invariant mass M of the diffractive system is above the prominent resonances which are excited 
diffractively from pions, specifically Ai(1260), 7r(1670), 7r(2100) and 7r(1300), 7r(1800) 
[111, ll^ and, perhaps, still higher radial and angular excitations of the pion. For instance. 



the color dipole model analysis in pO[ has shown that diffraction excitation of nucleons is 
exhausted by resonance excitation for M < 3 GeV. Therefore, the parton level calculation is 
viable at best for jets with fc > 1.5 GeV. 

The minor technical difference from diffractive excitation of the photon studied in [§, ^ 
|l^ is the change from the pointlike '~f*qq vertex ey4^\E'7^\E' to the non-pointlike nqq vertex 
ir7r(M^)^75\l'. In terms of the quark & antiquark helicities A the TTq{k)q{~k) vertex has the 
form (for the related discussion see Jaus [p]] ) 

^,(fc)75^A(-fc) = ^=L=[mj5,_-, - V2k ■ e.,6,-,] , (2) 

where ruf is the quark mass and ex = — (Ae^. + iey)/\^. In transitions of spin- zero pions 
into qq states with the sum of helicities A + A = ±1 the latter is compensated by the orbital 
angular momentum of quark and antiquark. 

In close analogy to the QCD description of diffractive dijet excitation in DIS, 7* -^ qq 
developed in 0, ||, |1^, the two helicity transitions in (|]) define the two diffractive amplitudes 
$o(^, k, A) and ^i{z, k, A). The lower blob in diagrams of Fig. 1 is related to the off-forward 
and off-diagonal differential gluon structure function of the target proton JF (xi, X2, k. A). In 
the considered high energy limit the two-gluon exchange interaction of qq states with the 
target conserves the quark and antiquark helicities exactly. This quark helicity conservation 
simplifies substantially the calculation of multiple Pomeron exchanges in diffraction off nuclei. 
One can readily update to the pion beam an analysis of the A dependence of diffractive 
amplitudes carried out for diffractive DIS in ||2^ , but for the purposes of diffraction on nuclei 
we only need the amplitudes for ttN -^ N{qq) in the forward limit A = and suppress A as 
an argument of diffractive amplitudes wherever it is appropriate. The lightcone momentum of 
the gluon is related to the change of the mass of the diffractive system, Xg = (M| — Mf^)/W^. 
In the considered problem Mf^^ = ml can be neglected, the diffractive mass M is generated 
by the first gluon, 

2;i~2;ip = — , (3) 

whereas the second exchange changes it only a little and X2 ~ 0. The detailed discussion of 
kinematics is found in 0, |^ and need not be repeated here, the principal point is that in the 
diffractive limit of Xjp ^ 1 the relevant off-diagonal differential gluon structure function of the 
target proton can be approximated H, by the conventional DGSF taken at x = |(xi +X2) = 
fa^ip, i.e., 

^ (xi, X2, /., A) = ^(ixip, /., A = 0) = ^^^f°^'f^ . (4) 

2 o log K^ 

After this rescaling one recovers precisely the expressions of ||^ for the diffractive amplitudes 
$2,^1- The hard scale in qq excitation is set by the large transverse momentum of jets. 



fc^ ^ 1 GeV^, and it is as{k'^) which enters the gluon-quark and gluon-antiquark vertices in 
the diffractive amphtudes. 

We find it convenient to introduce 



4^ /■ ,2,^(i3^iP'^^) 



and the distribution function 



f^'\k) 



fc^ 



47r J^Uxjp,k'^ 



SCTO 



fc^ 



(5) 



(6) 



normahzed to unity: j d?kf^^\k) = 1. For the sake of brevity of notations we suppress the 
dependence on xip. 

We define the two diffractive amphtudes ^o{z, k) and $i(-2, k) as (we use the normahzation 
shghtly different from that in [^, 0): 



$o(;z, fc) = as{k'^)(To I d^Kmf[ij^{z,k) -ij^{z,k- K)]f'^^\K) 

= as{k'^)(To d'^Kmfilj^{z,k)f'^^\K) — d'^KmfiljT,{z,K)f^^\k — K 
^iiz,k) = asik'^)aoJd^H.[k^^{z,k)-{k-H,)^lj^{z,k-K,)]f^^\K) 

= as{k^)ao fd^Kktlj^{z,k)f^\K,)-fd^KKiP^{z,K)f^\k-K 
The differential cross section of forward dijet production equals 



(7) 



da 



D 



dzdk'dA' 



A=o 



TV 
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{|<l>o|' + |$iP} 



(9) 



The radial wave function of the pion in momentum space is defined in terms of the nqq 
vertex function as 



ip^{z,k) 



(10) 



47r^z{l-z){M^ -ml) 

and is so normalized that the ir ^ fiu decay constant equals (we use the PDG convention 

F^ = 131 MeV ^) 

Fj, = d^kdzruf^^^z, k) = F^ dz(j)T,{z) . (11) 

Here we indicated also the relationship to the often discussed pion distribution amphtude 
(pTrlz) which for the purposes of our discussion we find it convenient to normalize to unity, 
J dz(f).^{z) = 1. We recall that for the pointlike photon r^(M^) = Cf, where Cf is the electric 
charge of the quark. In contrast to the pointlike photon for the non-pointlike pion r.,^{M'^) 
vanishes at large M^ faster than oc M~^, the relevant arguments are found in Brodsky & 
Lepage |]13[ and need not be repeated here. To this end we disagree with Ref. ^^ in which the 
pointlike r7r(M^) = const is assigned to the large-M^ tail of the pion wave function. 

Let us focus on the amplitude $o(^) k). The first term in (^ comes from diagrams la, lb, 
the corresponding spectrum of jets would be identical to the quark (antiquark) momentum 



distribution in the pion. Because the second term is a convolution of the gluon distribution 
and wave function, as such it is a broader function of k than V'ttI^, k) alone and would always 
take over at large k. The precise pattern of this dominance depends on the large- A; properties 
of ^(a^ip, fc^), the detailed discussion of which is found in [|T0|. Here we only mention that for 
X ~ 10~^ relevant to the E791 experiment the results of |1^ correspond to the inverse power 
asymptotics at large-A;^ 

/«(fc) oc A;-2^ (12) 

with the exponent 5 ~ 2.15. At s ~ 10~^ this asymptotics sets in at fc^ > A;^ ~ 1 GeV^, 
for smaller values of a^ip the exponent 5 is smaller, for instance 5 ~ 1.7 for xip = 10~^, and 
kf^ gets larger, see Fig. 2. As such, f^^\k) decreases at large k much slower than the pion 
wave function ip.„{z,k) (see the explicit parametrisation ( |iO| ) below) and the asymptotics of 
the convolution will be controlled by the asymptotics of f^^\k). We evaluate the second term 
in (!]) to the next-to-leading twist making use of the small-K expansion 



/«(fc-A.)^/«(fc) 1 



___5_ + 25(5 + 1)-^ 



(13) 



and obtain 

^o{z,k) =as{k'^)ao 
where 



TTli 



f^^{z,k)-f'\k)F^Mz)[l + 6- 



:{<{z)y 



e 



{<{^)) 



I d'^KK^-ipTr{z, k) 



(f'HlpT.^Z, n) 



(14) 



(15) 



In the related evaluation of the large-fc behavior of the convolution term in the diffractive 
amplitude $i for excitation of dijets with the sum of helicities A + A = ±1 the leading term 
comes from the second term oc 52kK,/k in the expansion (|T3l). Then we find 



^i{z,k) = as{k'^)aok 



^^(^,fc)_MWl/(i)(fc)F.0^(^) 
rufk 



(16) 



Evidently, in the region of large k where ipT^{z,k) dies out the amplitude #i will give the 
higher twist correction to the high-fc dijet cross section. 

The resulting large-fc asymptotics of the differential cross section for dijet production on 
nucleons reads 



da 



D 



dzdk'dA' 



A=o 



27r^ 



F:<Pliz)aUe) 



TUxjp^k) 



-\ 2 



fc^ 



26^ 



{<{z)) 



1 



2m) 



.(17) 



Evidently, the large-fc asymptotic behaviour oc k"^ is suggested by purely dimensional count- 
ing. Substantial departure from the law oc k~^ is possible because of scaling violations in 
^{\xjp,k). According to the recent phenomenological analysis |T^, the DGSF J-'{x,k) is 
approximately constant at moderately small x ~ 10~^, but rises steeply with k at x < 10^'^. 



Also, the experimental data are taken at fixed W'^, so that in view of (|l]) and (j^) the k^ depen- 
dence of the observed cross section is affected by the increase of xjp and decrease of J-'^^xjp, k) 
with increasing k . Similar kinematical bias affects the 2;-dependence of the experimentally 
observed cross section. 

There are three important aspects of our diffractive dijet excitation amplitudes at large k 
where the pion wave function dies out. 

First, here both helicity amplitudes are proportional to the DGSF of the target proton 
•^(|2;ip, A;^), i.e. the jet momentum comes from the momentum of gluons in the exchanged 
Pomeron, hence the term "splitting the Pomeron". To this end we recall that Nikolaev and 
Zakharov found the same proportionality of diffractive amplitudes to J^(|x]p, k ) also for real 
photoproduction with pointlike 'jqq QED vertex |^. From here one would conclude that 
this property does not require the wave function of the pion to be soft and the Trqq vertex 
function r^(M^) to vanish at large M^. Here we disagree with |]14[ who claimed that diffractive 
amplitudes are proportional to the integrated gluon structure function G{x,k'^). We notice, 
however, that in real photoproduction the cross section is dominated by the contribution from 
the helicity amplitude $i rather than $o in the pion case. Also, because of the pointlike 'yqq 
QED vertex the photoproduction cross section is oc k~^, see eq. (29) of H], compared to k~^ 



for pions as given by eq. ([T7|) . 

Second, in the same regime of hard dijets both diffractive amplitudes are proportional to the 
pion decay constant F^ and, more important, to the pion distribution amplitude (pni^). By the 
nature of our derivation this property emerges if the radial wave function of the pion ip{z, k) 
is a steeper function of k than f^^\k), which holds naturally for the anticipated decrease 
of non-pointlike irqq vertex function and for the phenomenologically known gluon structure 
function of the proton. Consequently, the ^-distribution of dijets allows the determination of 
the 2;-distribution of the pion distribution amplitude 4>n{z). 

Third, we emphasize that to the leading twist the differential cross section for dijet produc- 
tion on nucleons (eq.|l^ does not contain any free parameters, and thus is the perturbatively 
calculable quantity. 

3 Nuclear diffraction amplitudes 

We consider coherent diffraction 

7!- A -^ qqA' , 

where the recoil nucleus A' remains in the ground state. We focus on the forward diffraction 
cone A^ < R^^, where Ra is the nuclear radius. The longitudinal momentum transfer to the 
nucleus equals A^, = x-jprriN and coherent diffraction is possible if A^ -C Ra^, which condition 
is satisfied in the E791 kinematics in which xip ~ 10^^, see also the discussion in section 4.1. 
At xjp ~ 10^^ nuclear effects in DIS are dominated by nuclear shadowing of the qq Fock 



state of the photon p5|, ^, ^. Hence one must sum the qq multiple-scattering amplitudes 
of Fig. 3, we show a representative set for the impulse approximation, j = 1 (Figs. 3a,3b), 
and double scattering, j = 2 (Figs. 3c-3e). The typical multiplicity of rescatterings, j, is 



much smaller than the target mass number A. Because of the quark and antiquark helicity 
conservation one can calculate first the gg-nucleus scattering amplitude and convolute it with 
the pion wave function. Because the radius of nuclei Ra is much larger than the pion radius 
R-rr one can safely neglect the A dependence coming from the gg-nucleon scattering and take 
the gg-nucleon amplitudes in the forward limit A = 0. The strong coupling enters the qq 
loop as as{k'^)- In the high-energy limit of xjp <^ I/Raitij^j the calculation and summation of 
nuclear multiple scattering amplitudes is readily done in the impact parameter representation 
[|15| , |28| , ^ , ^ . Namely, we notice that after passing to the qq color dipole representation the 
helicity amplitudes ^o{z, k) and ^i{z, k) can be cast in the form 

j2„ „—ikr 



%{z,k) 
^i{z,k) = - 



dre a{x,r)mf-^^{r,z) 

(fr e~^^^ a{x, r) V^^(r, z) . 



Here k is the transverse momentum of the jet, r is the qq separation in the impact parameter 
plane, 

'^Jz, r) = t::^^^ I (fkip^iz, k) exp{ikr) (19) 



is proportional to the qq color dipole distribution amplitude in the pion and 

a{x, r) = as{k'^)(To f (fK,f^^^ (k) [1 - exp(iKr)] (20) 

has the meaning of the dipole cross section for interaction of the qq dipole r with the target 
nucleon in which the strong coupling as enters at the hard scale given by the jet transverse 
momentum k. 

The Glauber-Gribov representation of amplitudes (|18|) for the nuclear target is obtained 
by substitution of the gg-nucleon scattering amplitude by the gg-nucleus scattering amplitude 



I BHI and reads: 



$i^)(^,fc,A) 



2m i 



2m I 



-ibA — ikr 



'^n{z,r) h -exp 



^a(x,r)T^(b) 



d'b / d're- 



vl/,(z,r)^(-i; 



n+l 



n>l 



a [x,rj 
2'^n! 



T2{b) 



^^{'\z,k,A) 



-2i / d'b / d're 



-ibA — ikr 



E(-ir 



+1 



a"(a;,rl 



n>l 



2'^n! 



nib) 



V^,(z,r), 

(21) 

where b is the pion-nucleus impact parameter, T^ib) = J dz'nA{b,z') is the familiar nuclear 
optical thickness JlSl, nA{b, z') is the nuclear matter density. The frozen dipole approximation 
(^TD is applicable because A^ <^ R'^ (HH, for the modern formalism see |2^). For the sake 
of simplicity above we took the exponentiated form for the nuclear profile function instead of 
its more exact form 



^A{b) 



a{x,r)TA ib) 
2A 



-\A 



1 — exp 



--a{x,r)TA{b) 



(22) 



reformulation of all results for the polynomial form poses no problems. 

Representative diagrams for the impulse approximation, j = 1, are shown in Figs. 3a, 3b. 
They give the familiar result 



<^i''\z,k,A) = <l>o{z,k) J d'be-'^^TA{b) = A<l>o{z,k)GeU^) 



(23) 



where Gem{A) is the charge form factor of the nucleus. 

The principal effect of rescattering is readily seen from the double scattering, j = 2. Making 
use of the integral representation ( ^Of ) we find 



/ 



d^re ^^'^^^{r,z)a^{x,r) 



a's{k')at, d'K,d'K,f'^{K^)f'^{K, 



f(i), 



d r e 



-ikr 



= al{k')al j d^,^,d^,^,f'\,^,)f'\,^,) 

X [^^(fc, z)-2 m^{k - ni,z) + ^^(fc - Ki - K2, z)] 

= aUe)al [vl/,(fc, ^) - 2 1 d-'K^^iH,, z)f^^ (k - h,) 



'^Ar,z) 



where 



+ j d'^n^^{i^,z)f^\k-n) 



f^\k) = Jd'K,d'f^2f'\K^)f'\K2)6{k - /^i - K2) 



(24) 



(25) 



is normalized to unity: J d'^kf^'^^k) = 1. The sum of the impulse approximation, n = 1, and 
double-scattering terms equals (for the sake of illustration we take A = 0): 



1 - -as{k^)aoTA{b) 



<l>i^\z,k,A = 0) = as{k')a^mf j d^hTA{h) [ip.{z,k) 

- f d'^Kllj^{z,K)f'^^\k-K) 

- ^ J d^KMz,i^)f^^\k - K)as{k^)aoTA{b)^ . 



1 - -as{k^)(ToTA{b) 



(26) 



The three terms in the last line of ( P4D and in (^) correspond to the three classes of double 
scattering diagrams shown in Figs. 3c-3e. The first term in the r.h.s. of (^6]) shows that 
the no-Pomeron splitting term in (|1^) coming from Fig. 3a in the nuclear case receives the 
conventional shadowing correction from the double-scattering diagram of Fig. 3c (and the 
not shown here partner diagram in which the two gluons from second nucleon couple to the 
antiquark). The Pomeron splitting term in (|T^) coming from Fig. 3b in the nuclear case is 
similarly shadowed by double-scattering diagrams of Fig. 3d. The effective shadowing cross 
section equals 

ae//(fc') = as{k'')ao (27) 

The new feature of double scattering is the third term in ( P^D given by the double-Pomeron 
splitting diagram of Fig. 3e. The convolution (|25| ) implies the broadening of f^'^\k) compared 



9 



to f^^\k). Furthermore, this broadened distribution /(^)(fc) has the same sign as, i.e. it is 
an antishadowing correction to, the single-Pomeron sphtting term (|26|) and would eventually 
take over single-Pomeron splitting at large k. 

Higher order rescatterings give rise to distributions with j-fold Pomeron sphtting 



•' i=l 



(28) 



which obviously broaden with increasing j. Rearranging nuclear diffractive amplitudes as an 
expansion over f^^\k) we obtain 



<''\^,k) 



2m/ Y. / d^^^' 



X 



a,jfik')TA{b) 



J 
2mf 



exp 



TAih) 



^effik'^) 



d^be-'^^\^^{z,k) 



2 
1 — exp 



^effik'] 



TA{h) 



Y^fd'K^^{z,K)f^^\k-K)^ 
j-yl-l J- 



ae//(fc')T^(b) 



exp 



0", 



eff 



'^ 



-TA{b) 



^['^\z,k) = 2j2fd^be-''^^fd^K[k^^{z,k)-{k-K,)^^{z,k-K,)]f''^\K) 



X 



a,ff{k'')TA{b) 



u 



exp 



cr. 



eff 



k' 



d^be-'^^ 



k'^^{z,k) 



2 
1 — exp 



TAib) 



J2 fd^^K,^il^iz,K)f^^\k-H,)^ 



^T.(6)) 

a,ff{k')TA{b)y 

exp 



(^eff{k^] 



■TAib) 



(29) 



The nuclear attenuation factors show that shadowing is indeed controlled by aefjik^)- Despite 
the decrease aefjik"^) oc ^^(fc^) numerically this cross section is quite large, grows slowly at 
very small xip, and is a soft gluon exchange dominated quantity. 

At large jet momentum the diffractive amplitude is dominated by the second term in (|29| ) 
in which all broadened j-Pomeron splitting contributions enter remarkably with the same 
antishadowing sign. 

Whether this antishadowing takes over shadowing depends on the large k behaviour of 
f^'^\k). For the power asymptotics (|12D the leading contribution to convolution f^'^\k) (p5[) at 
large k comes from the configurations when there is one hard splitting of the Pomeron with 
Ki r^ k, whereas all other ki are small. Correspondingly, the values of gluon lightcone momenta 
Xi in the hard splitting of the Pomeron are the same as in the free nucleon case, whereas in the 
predominantly soft rescatterings x, ~ X5 ~ IGeV'^/W'^ and in multiple scattering expansion 



(|29D all aeff{k ) but one for the hard splitting of the Pomeron must arguably be evaluated at 



10 



X = xs- However, we notice that according to |T^ the x-dependence of J^{x, fc^) is weak for 
soft k^ and in the practical evaluation of f'^^\k) one can put xs = x^l2. Then to the leading 
twist one readily finds the large-fc^ asymptotics 



PKk)=jf^'\k) 



Wi 



(30) 



which clearly shows the anticipated enhancement of the large-fc^ tail by multiple rescatterings. 
The salient feature of large-fc broadening (pOD is that the exponent of the power asymptotics 
is sustained. Because of the normalization condition j d?kf'^^\k) = 1 the broadening at large 
k entails the small-fc depletion f^^-'^\k) < f^^\k) at small-fc. Evidently, for larger j this 
small-fc depletion will be stronger and would extend to larger k. 

The evaluation of the higher twist correction making use of the expansion ([l^) proceeds as 
follows. For the sake of definiteness focus on the contribution from the configuration in which 
Ki ^ k and the total transverse momentum for the subset i = 2, ..., j is small, {J2i>2 f^if ^ ^^• 
Then, to the Leading Logfc^ accuracy the higher twist contribution will be dominated by the 
(j — 1) configurations of the subset 2 = 2, ..., j in which one of the nf is running up to ~ k^ 
whereas all other momenta are small: 



fU){k)c:.jf^'\k)Jd'K2...d'^, 



k' 



l + 7^(E'^i 



Jf^'\k) 



1 + ^(J-1)/ d^'^^'f^'^ 



K 



3r\k) 



1 + 



i>2 



/«(a.2).../«(a^,) 



(j-l)G'(-xip,fc2 



ScTofc^ "^ ' ^2 



(31) 



Remarkably, the coefficient of the higher twist correction is proportional to the gluon structure 
function of the proton G{^xjp^k^) at a hard scale k^. For the j-Pomeron splitting diagrams 
the higher twist correction rises oc (j — 1), i.e. the antishadowing/broadening of large-fc^ tail 
takes place also to the higher twist. 

Making use of the result (|3TD in the expansion (29) we obtain 



■M) 



$^^^(2,fc,A) = 2 / d'be 



-ibA 



ruf^^ik, z) 1 - exp ( --CTe//(fc')TA(6) 



--/«(fc)F^0^(^)aoa5(fc') 



T.(6) + ^^i^G(.^)Tl(6) 



(32) 



Here the impulse approximation term oc ^T^ik^z) is shadowed with the soft cross section 
(pTl). The most remarkable feature of the Pomeron splitting contributions (^) is an exact 
cancelation of soft shadowing and antishadowing/broadening effects. Furthermore, this can- 
celation makes redundant the exact value of xs at which ctq must be taken in the nuclear 
multiple scattering expansion. The broadening law (^) is crucial for this cancelation of soft 
shadowing and antishadowing/broadening effects. Similar exact cancelation of shadowing and 
antishadowing/broadening effects, and independence on the exact value of soft shadowing cross 
section o"o, take place in the next-to- leading twist correction too. Consequently, both leading 
and next-to-leading twist amplitudes are parameter free calculable in hard perturbative QCD. 
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Following the derivation (IC) for the free nucleon target, a similar analysis can be repeated 
for the nuclear target with the result 

1 



.(^) 



^r{z,k,A) 



2k 



d^be-'^^\^^ik,z) 



1-exp --ore//(fc^)T^(6) 



- l:aoas{k')^-^^^i^f^'\k)F^M^)TA{b) 

I VTifk 



f33) 



where we neglect the corrections oc /;; ^ to the already higher twist convolution term. 
We notice that within the diffraction cone of A^ < i?^ the viable approximation is 



She-'^^Tlih) 






(34) 



where the coefficient Ca ~ 1 depends slightly on the shape of the nuclear matter distribution: 
Ca = 1 for the Gaussian density appropriate for light nuclei decreases slowly to Ca = ^ for 
the uniform density sphere. Then for large jet momentum k when the impulse approximation 
contribution dies out, our nuclear diffractive amplitudes (0), ( ^3]) take a particularly simple 
form 

Anas{k')J'{\xTP,k'') 



$S,^)(z,fc,A) = -AF^0^(^) 



k" 



1 + 6 



2^^li^)y 



e 



X G., 



emK^l + 



2{Rl,)e G{---'^k)GU-,^) 



(35) 



#S^)(^,fc,A) 



-kAF^cP^iz) 



4Tras{k'^)J^{lx-sp,k'^ 
3 F 



"Cemf A 



2^5KM 
nifk 



(36) 



Finally, a correction for the finite longitudinal momentum transfer to the nucleus can be 
evaluated as follows. First, within the diffraction cone the impulse approximation ampli- 
tude ( p3D acquires the longitudinal form factor Gem(3^n'"^Af)- Second, we have shown that 
in multiple rescattering contributions there is only one hard splitting of the Pomeron with 
the longitudinal momentum transfer A^ ^ xjpmN whereas in soft rescatterings the longitu- 
dinal momentum transfer can be neglected. As a result, the same longitudinal form factor 
G emix^m^j) holds for all multiple rescatterings. 



Upon the A integration we find that the large-fc asymptotics of nuclear diffraction cross 



section will read 

dan 



dzdk 



27r5 



Fyjz)GUxl.m%)aUk' 



•^(o-^IP) "^ 



2An2 



k' 



3A^ 



2{^liz)) 1 



(k2(^))\ , 2nGAAas{k' 



2m) 



^RIh) 



-G{-xjp,k'^ 



(37) 



where we included the effect of the longitudinal form factor. Clearly, at a sufficiently large k 
the higher twist correction will be dominated by the last term which is model-independent, 
enhanced by the gluon structure function G{hx-sp, k ) and rises for heavy nuclei oc A^. 
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4 The numerical results and comparison with the E791 
data 

4.1 Nuclear broadening of multiple-Pomeron splitting contribu- 
tions 

In the E791 kinematics 1.25 GeV< k < 2.5 GeV and W"^ = 940 GeV^, so that xip ~(1-2)-10-2. 
We start with a check of the accuracy of the expansion (|3lD according to which 



must exhibit j-independence at sufficiently large k . The results for xjp = 2 ■ 10~^ shown in 
Fig. 4a demonstrate this is indeed the case for very large k > (2-3) GeV. Slight departure from 
the universality can be understood in terms of the slight k dependence and the broadening 
driven j dependence of the effective exponent S. The situation at lower fc is a nontrivial one, 
because at a sufficiently small k the broadening of f^^\k) must be superseded by the small-fc 
depletion, in which region 

Aj{k)r~.--. (39) 

The effect of depletion extends to larger k with increasing j. Indeed, for k = 1 GeV and 
j = 10 our Aj{k) is getting close to the no-broadening estimate (^). Anyway, the finding of 
Aj(/c) < indicates breaking of the broadening law (|30| ) for j > 5 at A; ~ 1 GeV, the point of 
crossover Aj{k) = is moving to larger k with increasing j. The large-/c asymptopia is even 
more elusive at smaller x, see Fig. 4b for xjp = 2 ■ 10~^. Here the broadening law ( |5D| ) is only 
applicable at A; > (4-5) GeV. 

According to these results, the E791 range of k falls in the transient region in which 
the higher twist expansion (|3lD is not applicable yet and the multiple-scattering broadening 
is weaker than given by (0). This suggests that in this transient region of k the large- /c 
broadening of f^^\k) is not sufficient for exact cancelation of shadowing effects. All numerical 
estimates of the A-dependence must use exact numerical results for f^^\k). 

4.2 The pion wave function and z-distribution amplitude 

In numerical calculations of diffraction amplitudes we use a slight modification of the Jaus 



2^] parameterization of the pion wave function, which in our convention (p!0| ) is 

1 / 1 



i^.iz.k) (X ^^^ _;^^, exp [--Rl{M' - Am))) (40) 

and with i?,^ = 2.2 GeV~^ and nif = rriu^d = 0.215 GeV provides a consistent description 
of the IT ^ fiu decay constant F.,^, charge radius of the pion, 7r° -^ 27 decay rate and slope 
of the 71^ -^ 77* form factor [^ (for the related analysis see [^). The numerical results 
for (j)Tr{z) are shown in Fig. 5 and differ only weakly from the often discussed asymptotic 
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distribution amplitude (j)asym{z) = Qz{l — z) (for the review see [0, |T^, |^), in the broad 
range 0.2 <z < 0.8 the difference does not exceed 10%. A convenient analytic approximation 
to (f>TT{z) given by this soft wave function is 

,,,.). 0,6572log(l + ML=i))e.p(-^) (41) 

and is good to better than 1 per cent apart from z < 0.03 and 1 — z < 0.03. 

The large-fc asymptotics of the hehcity amplitude for excitation of dijets with A + A = ±1 
is proportional to (k^(z)), defined in ([T5|) . The soft wave function (^) gives 

i-liz)) ^ P^^ ^ , (42) 

R log 1 + p2 2 

i.e., (k^(-2 = 0.5)) ~ 0.17 GeV^, which is a natural scale for the soft pion wave function. It 
decreases gradually away from z = 0.5, for instance {kI{z = 0.2)) = {^^{z = 0.8)) ~ 0.12 
GeV^. 

Here we notice that in view of (|1|) and (^ there is the kinematical z-xjp correlation xip oc 

[4:z{l — z)]"^. The xip dependence of DGSF can be parameterized as J-'^^xjp, fc^) oc Xjp . 
With allowance for the z-x^p correlation the observed 2;-dependence of diffractive amplitudes 
changes from 07r(-2) to 

0f(r,^)(X0.(z)[^(l-^)r(^'). (43) 



The phenomenological determination of the exponent r(fc ) in [|T0| gave r ~ 0.16 at /c = 1.35 



2> 



GeV and r ^ 0.22 at k = 2 GeV for the GRV-D and MRS-D parameterizations for JF(ixip, k 
and r ~ 0.25 at k = 1.35 GeV and r ~ 0.30 at k = 2 GeV for the CTEQ-D parameterization. 
In Fig. 5 we show the observed z-distribution amplitudes for r = 0.15 and r = 0.30, evidently 
the observed z distribution gets even closer to 4>asym{z). 

4.3 Importance of multiple Pomeron splitting processes in nuclear 
diffraction 



The result for the Pomeron splitting term in ( ^2] ) has a deceptively simple form of the sum 
of the single and double scattering terms but such an interpretation would be utterly wrong. 
An importance of multiple nuclear rescatterings can be judged from the mean value of j in 
expansion (|25|). Casting it in the form $g '{z, k) = J2j uijf^^\k), for k =2 GeV we find 

,. ^ E.w.jf^Hk) ^ r 2.03 for ^'C 

EjWjf(^\k) \ 4.10 fori96pt^ ^ ' 

which shows clearly an inadequacy of truncation of nuclear rescatterings to the single and 
double scattering. Indeed, at x ~ 10"^ and fc ~ 4 GeV^ the effective shadowing cross section 
is quite large, a^/fik'^) ~ 40 mb. Closer inspection shows that for exhausting 95 % of the 
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strength of $^ one needs the contributions up to j = 4 for the carbon target (A=12) and 
up to j =(8-9) for the platinum target (A=196). For the E791 energy, i.e. xjp ~ 10~^, this 
imphes in conjunction with the results shown in Fig. 4 that from the viewpoint of convergence 
of nuclear expansions k > 2.5-3 GeV are needed for the applicabihty of the leading plus 
next-to- leading expansions (p5D, (|36D, (pT]). 



4.4 Evaluation of next-to-leading twist corrections 

The next-to-leading twist correction can conveniently by parameterized as H/k"^. First we 
evaluate the contribution to H coming from the pion wave function. It is controlled by the 
moment {kI{z)). For the 0.2 < z < 0.8 which are relevant to the E791 data, we take in further 
estimates {k,'^{z)) = 0.15 GeV^. We notice that numerically {k'^{z = 0.5)) ~ 3.5m'j, which 
entails that on the free nucleon target next-to-leading twist effects come predominantly from 
the helicity component oc $^ in the cross section, see eq. ([T7|): 

hW = fM(^ ^ 0.555^ GeV-2 ~ 2.2 GeV-' . (45) 

mj 

Here we took 5^ ~ 4 appropriate for k > (2.5-3) GeV where the asymptotic expansion (|37D is 
applicable. 

A similar estimate for the contribution to H from from the helicity component oc $0 i^ 
the cross section is 

i7|°) ~1.3GeV"^ (46) 

For xjp ~ 10~^ and k' ~ (3-4) GeV^ relevant to the E791 kinematics, the analysis [|10| 
gives asG{^x-sp, fc^) ^ 1-1.2. Then the nuclear rescattering/broadening contribution to next- 
to-leading twist is (we use the nuclear density parameters from the compilation p9[| ) 



AtiS'CaA' „2^^/l ,2N /n.. noN.2.1 fl.5GeV2 for ^^C 
^^^^«.(/.^)G(-x., e) . (0.16 - 0.2)S'A^ ^ 3_^ ^^^, ^^^ ,,,^^^ 



Lumping together all three contributions, we find H^^ ~ 7 GeV^ and if" ~ 5 GeV^. Conse- 
quently, no simple expansion in the leading and next-to-leading contributions is possible in the 
E791 region of jet momentum and an accurate numerical evaluation of broadening of f^^\k) 
is called upon. 

4.5 Comparison with the E791 data: k distributions 

In Fig. 6 we show our numerical results for the A^-integrated nuclear diffraction cross section 
for the platinum target. We are interested in the large k where the contributions oc ip^iz, k) 
died out, which in both helicity amplitudes is preceded by the zero. In order to have a crude 
idea on where this happens, we stretched our calculations for platinum target down to k = 0.5 
GeV. The amplitude $0 has a zero at A; ~ 0.65 GeV. Because the term oc ipTri^^ k) is sensitive 
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to nuclear shadowing and to the soft cross section ctq thereof, the position of the zero is model- 
dependent. In $1 the second component is of higher twist and $i has a zero at larger fc ~ 1 
GeV. The impact of these zeros of diffractive amplitudes is manifest up to k > 1.5 GeV, 
which is still another reason why a comparison of our predictions with the experimental data 
is justified only at k > 1.5 GeV. 

The E791 data give only the k dependence of the acceptance-corrected cross section without 
absolute normalization. The normalization of our theoretical curve is the eyeball fit to the 
data, the agreement with the experimentally observed ^-dependence is good at fc > 1.5 GeV. In 
the theoretical calculations we include consistently the k -xjp correlation discussed in section 
2 following eq. (|17D . 

At E791 energy our leading twist theoretical cross section would have followed the law /c~" 
with the slope 

n ^ 4 + 2(5 + 2r(fc2) ^8.7-8.8. (48) 

Our numerical results do not exhibit simple dependence ex fc~", if we define the local slope n 

as 

d log an 

dlogk^ 

then for k > 2.5 GeV we find n ~ 12, whereas around k = 2 GeV the slope n ^ 10. We 

attribute these large values of local n to very large higher twist corrections. 

In Fig. 6 and Fig. 7 we show separately the contributions to the diffractive cross section 
of the two helicity states: the leading plus next-to-leading twist component with A + A = 
and the higher twist component A + A = ±1. They are of comparable magnitude, the higher 
twist component oc $^ starts dying out and the component oc $q starts taking over only at 
k > 2.5-3 GeV, in perfect agreement with evaluation of higher twists in section 4.4. 

Fig. 7 shows our predictions for the future experimental tests of the absolute normalization 
of diffractive cross section. In all the numerical calculations we used the parametrization 
labeled 'D-GRV from [l^. Here for the sake of convenience we plot k^ da d / dzdk'^\z=Q.^- As an 



illustration of the energy dependence of dijet production here we also show the predictions for 
Ej, = b TeV. In this case xtp ~ 10~^ and 5 ~ 1.7, so that in the leading twist we expect n ~ 8. 
Furthermore, because of the lower value of 5 the higher twist effects would be about half of 
those for E^^ = 500 GeV in the E791 experiment, cf. Fig. 4a and 4b. Indeed, for E^ = 5 TeV 
k^da£)/dzdk'^\z=o.5 exhibits much weaker dependence on k^. 

4.6 Comparison with the E791 data: nuclear mass number depen- 
dence 

The A-dependence of nuclear cross sections is often parameterized as a oc A". The E791 
experiment uses the carbon and platinum targets and defines the exponent a as 

log — 

a = ^ 



1 ^1 
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For the reference, in the impulse approximation 



a oc 



(^c\) 



and simple evaluation for Pt and C nuclei gives ajA = 1-44 (the average slope for the A 
dependence in the same interval gives ajA = 1.39, for instance, see |3^). A comparison of 



our numerical results for Pt and C nuclei yields the /c-dependence of the exponent a shown in 
Fig. 8. 

As we discussed in section 4.1, at k < 1.5 GeV antishadowing/broadening effects are too 
weak to cancel nuclear shadowing, which explains the small value of a{k ) < ajA- On the 
other hand, at larger k higher-twist corrections become substantial and give a{k'^) > (Xia at 
k > 1.8 GeV. At large h we find again good agreement with the E791 results. 

4.7 Comparison with the E791 data: ^-distributions and the pion 
wave function 

To leading twist diffraction of pions into dijets uniquely allows to measure the pion distribution 
amplitude 0,1 (-z). As we have seen above, for moderately large k studied experimentally higher 
twist corrections are very large. The parameter {k'^{z)) of higher twist correction from the 
pion wave function varies with z and is a model dependent one, our estimate (|^) gives only a 
crude idea on its z dependence and (|^), (|46|) must be regarded as numerical estimates within 
the factor two. Notwithstanding these uncertainties, even on the free nucleon target one must 



be able to distinguish experimentally between the double-humped CZ |TT| and asymptotic 
distribution functions. 

To the contrary, the nuclear-rescattering driven higher twist correction is model- independent. 
The xjp-z correlation driven z dependence of G{^x-sp, fc^) is very weak, it only slightly enhances 
the cross section around z = ^. For this reason, even if higher-twist dominated, the diffraction 
off heavy targets is a good probe of the pion distribution amplitude (pTr{z). The E791 paper 
does not give the acceptance-corrected 2;-distributions. Ashery concludes 0] that within the 
~ 20% experimental error bars the observed E791 ;2-distribution is consistent with the Monte 
Carlo modeling based on the asymptotic pion distribution (pasymiz). Because the observed 
^-distribution (^) given by our model soft wave function is very close to the asymptotic one, 
it is perfectly consistent with the E791 data. 

5 Conclusions 

We developed the perturbative QCD description of diffraction dissociation of pions into hard 
dijets on nucleons and nuclei. To leading twist dijet excitation is shown to be dominated by 
the Pomeron splitting mechanism and the two diffractive helicity amplitudes are shown to 
be proportional to the unintegrated gluon structure function of the proton. We derived an 
multiple-Pomeron splitting expansion of nuclear amplitudes which is of antishadowing nature. 
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To leading twist there is a remarkable cancelation of nuclear attenuation and antishadow- 
ing/broadening effects. We obtained a model-independent estimate for next-to-leading twist 
corrections driven by nuclear rescatterings. These higher twist corrections are shown to be 
very large up to jet momenta k > (5-7) GeV^ and affect substantially the k and atomic 
mass number of the diffraction cross section. The model dependence of extraction of the z 
dependence of the pion distribution amplitude is shown to be weak for diffraction off heavy 
nuclei. Our calculations based on the recent determination of unintegrated gluon structure 
function of the proton reproduce well the basic experimental findings from the recent E791 
experiment. A simple interpretation of the observed k dependence is not possible though 
because in the E791 range of k the cross section is overwhelmed by higher twist effects. 
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Figure 1: Feynman diagrams for diff'ractive dijet excitation in ttN collisions. 
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Figure 2: The k-distribution f'^^\k) as a function of k for several values of xjp 
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Figure 3: The nuclear multiple scattering series for diff'ractive dijet excitation on nuclei. Di- 
agrams a),b) are sample diagrams of impulse approximation, diagrams c)-e) represent the 
various types of double scattering contributions. Higher order contributions, that appear in the 
calculation are not shown. 
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Figure 4: Large k-scaling properties of the multiple convolution integrals. Shown is Aj{k) for 
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Figure 5: The pion distribution amplitude (p-K^z). The solid line is (p-^^z) calculated from the soft 
wave function eg. (|7^j, the dashed line is the asymptotic distribution amplitude. The curves 
labeled 'observed' show the z-dependence of the soft pion distribution amplitude modulated 
with the z-dependence due to the kinematical xjp — z correlation for two different values of 
the effective exponent t, see eq. (Oj 
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Figure 6: The E791-data ^ for the differential diffractive dijet cross section da/dk for the 
^^^Pt target with the theoretical calculations. The data are not normalized. The dash-dotted 
line shows the contribution of the helicity amplitude $o ; ^he dashed line is the contribution 
from ^i . The solid line is the total result. 
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Figure 7: Theoretical predictions for the differential dijet-cross section for the ^^^Pt target. 
Panel a) is for the energy of the E791 experiment, E = 500 GeV , panel b) for E = 5TeV. The 
dash-dotted line shows the contribution of the helicity amplitude $0 , the dashed line is the 
contribution from ^i . The solid line is the total result. 
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Figure 8: Exponent a of the atomic mass number dependence of the dijet cross section with 
the results from E791 j^. The dashed line shows the impulse approximation result, the solid 
line is the result of the full multiple scattering series. 
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